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Abstract: 


The growing interest in Laser Additive Manufacturing (LAM) High-entropy alloys (HEAs) during 
most recent years, the design of elements and scanning strategies are primary methods to overcome the 
defects such as porosity, pores, and segregation etc. Here we propose a new approach, a novel real-time 
Laser Shocking of Melt Pool (LSMP), to obtain melt pool modifications to yield better HEAs. LSMP 
moves a pulsed laser into a liquid melt pool caused by a continuous wave laser, enabling 
non-destructive and real-time fabrication of high-performance HEAs. The numerical simulation reveals 
the convection mechanism of the melt pool in the LSMP process, and the intervention of the pulsed 
laser promotes melt pool flow type to convert the Marangoni effect into a multi-convective ring, which 
accelerates melt pool flow and the columnar crystal growth is significantly inhibited. Experimental 
results show the evolution law of the microstructure in the LSMP process. CrFeCoNi HEAs forms a 
stable FCC structure without segregation and other defects. Meanwhile, the microstructure undergoes a 
Columnar-Equiaxed Transition (CET), and CrFeCoNi HEAs with higher nano-hardness are obtained. 
Furthermore, the reduction in the (111) interplanar spacing demonstrates laser shock induces beneficial 
compressive stress. Therefore, it is necessary to consider introducing laser shock at appropriate 


conditions to control the melt pool online. 
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1. Introduction 


As a new type of material, Yeh et al. [1] broke the traditional design concept of a single metal as 


the main element and created High-Entropy Alloys (HEAs) with excellent performance. It has brought 


a new dawn for the development of new alloy systems, which have received extensive attention in 


many fields of basic research and engineering manufacturing. But, typical metallurgical defects such as 


porosity, pores, and segregation in conventional as-cast HEAs affect the final properties, which limit 


the potential application of HEAs. Fortunately, Brief et al.[2] firstly used Laser Additive Manufacturing 


(LAM) technology to obtain HEAs with better mechanical properties than as-cast samples, 


demonstrating the feasibility of LAMed HEAs. It not only shortens the process cycle but also provides 


a faster and more efficient method for manufacturing complex HEAs parts. Besides, Moghaddam et 


al.[3] presented a conclusion of current state and future prospects of additive manufacturing (AM) of 


HEAs. 


Among many HEAs systems, laser additive manufacturing FeCoCrNi-based HEAs is the earliest 


and most extensive. Zhou et al. [4] fabricated carbon-containing FeCoCrNiCo9; HEAs by Selective 


Laser Melting (SLM). The alloy structure is a single face-centered cubic (FCC) without carbide phases, 


and the carbon distribution in the samples prepared under optimal parameters is uniform. Zhu et al.[5] 


prepared the CoCrFeNiMn HEAs by SLM, the distribution of each element is uniform, and Mn is 


slightly segregated at the melt pool boundary. Besides, Peyrouzet et al.[6] used SLM process to obtain 


Alo3;CoCrFeNi HEAs with a tensile strength of up to 896 MPa and a yield strength of 730 MPa, which 


is also a single FCC structure. Obviously, more and more scholars such as Li et al. [7] have transferred 


the research focus on HEAs to the performance of materials. 


However, the high thermal stress caused by the thermal cycle of the melt pool during the LAM 


process also brings many hazards. Zhang et al.[8]observed cracks in the CoCrFeNiMn HEAs samples 


obtained under different scanning strategies, resulting in low densities. Niu et al.[9] also found the 


same result. And many scholars have joined the ranks to solve this problem. For example, Jiang et al. 


[10] shows the presence of FCC phase significantly reduces the internal stress of the specimen by the 


addition of Fe-based amorphous. In addition, the indeterminate cooling rate of the melt pool, a critical 


factor, also directly affects the recrystallization species, leading to various types of HEAs with different 


properties. Li et al.[11] reveals the formation mechanism of microstructures in HEAs rapidly cooled to 


room temperature at nanoscale through Molecular dynamics (MD) simulation. At present, most HEAs 


prepared by Laser Melting Deposition (LMD) exhibit a typical dendritic structure, because the scan 


rate of LMD is significantly lower than that of the SLM process, causing a low cooling rate so that 


equiaxed crystals cannot be achieved. Therefore, Tong et al.[12] found the microstructure of the 


CoCrFeNiMn HEAs fabricated by LAM is mainly columnar crystals inside the melt pool, while the 


outside of the melt pool is equiaxed grains. After heat treatment at 1100 °C, the structure is completely 


equiaxed grains. Certainly, subsequent heat treatments such as Knowles et al.[13] used a range of 


annealing processes after Direct Metal Laser Sintering (DMLS) and Li et al.[7] applied the Hot 


isostatic pressing (HIP) after SLMed HEAs are detrimental to the continuity and production efficiency 


of the AM process. However, Tong et al.[14] shows different process parameters such as laser power 


and scanning speed in the LAM process can directly obtain different microstructures because the ratio 


of the thermal gradient and solidification rate inside melt pool is the dominant factor in 


Columnar-Equiaxed Transition (CET). 


Thence, the online regulation of the melt pool to get the better structure in the LAM process is 


particularly important. The Laser Shock Peening (LSP) as a advanced technology which not only 


induced high compressive stress on the surface of metal components and structures but also achieved to 


control at the micro-nano level. For example, Hu et al.[15, 16] realized ultrafast direct fabrication of 


plasmonic nanoarrays by LSP. A novel synchronous Laser Shocking of Melt Pool (LSMP) method is 


proposed to suppress thermal stress, enhance melt pool convections and inhibit columnar crystal 


growth. Not only that, it extends the effort that Tong et al.[14] used the LSP as a post-processing 


treatment to in-situ modification to improve the LAMed HEAs. Although the pulsed laser is widely 


used in material properties strengthening, to our knowledge there had not been any research on 


synchronous hybrid additive manufacturing processes with laser shocking of melt pool to overcome the 


surface tension of liquid metal.. In this paper, the action mechanism of LSMP by combining 


experiments and numerical simulations is studied. We reveal the convection mechanism inside the melt 


pool with the assistance of pulsed laser. As a result, the microstructure evolution law of HEAs 


produced by LSMP is obtained, achieving non-destructive real-time production of higher-performance 


CoCrFeNi HEAs. 


2. Material and Methods 


2.1. Materials and Experiment 


The LSMP method moves the pulsed laser into a melt pool of liquid high-entropy alloy produced 


by a continuous wave (CW) laser, as shown in Fig. 1(a). A size of 30 mm x 30 mm x | mm 45 steel, 


cleaned with ethanol solution after polishing to ensure cleanliness, was used as the substrate. The 


experimental material adopts the equiatomic CrFeCoNi powder prepared by vacuum atomization. The 


particle size is 15~45um. It has good fluidity and sphericity. The powder morphology and composition 


analysis are shown in Fig. 1(b) and (c). The thickness is 1mm. The experimental equipment includes a 


1kW continuous wave (CW) fiber laser (Shenzhen Chuangxin Laser Co., Ltd.) and a pulsed laser. The 


Selective Laser Melting (SLM) process is carried out in a high-purity argon atmosphere with a gas 
supply rate of 3L/min. The process parameters are laser power 800W, platform moving speed 
1000r/min, and spot diameter 3 mm. The pulsed laser process parameters are laser energy 3J, pulse 
width 10ns, repetition frequency 10Hz, and spot diameter 1mm. The CW laser beam is perpendicular to 


the substrate, and the pulsed laser beam deviates from the vertical by -10°. 
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Fig.1. (a) Experimental schematic diagram, (b) the SEM image of CrFeCoNi powder, (c) the EDS analysis of CrFeCoNi powder. 


2.2. Test Methods 


Different samples were tested from three perspectives: molding process, microstructure and phase 


composition to mechanical properties. First of all, the Chronos1.4 high-speed camera (Kron 


Technologies, Canada) was used to monitor the manufacturing process, and the NewViewTM 9000 


Series white light interference 3D profiler (ZYGO, The United States) was used to measure the profile 


of the sample. Secondly, the sample was moved along the normal direction, then it was corroded by 


aqua regia (hydrochloric acid: nitric acid=3:1) for 60 seconds after polishing and cleaning the 


cross-section. The microstructure and the Energy-dispersive Spectroscopy (EDS) tests passed the BX51 


optical microscope (OLYMPUS, Japan) and the MIRA3 scan electron microscope (SEM) (Tescan, The 


Czech Republic). To understand the phase composition of the original CoCrFeNi powder and the 
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samples, the D8 ADVANCE X-ray diffractometer (XRD) (Bruker, The German) was used for detection. 
The target material was copper and the scanning angle was 20-90°, the scanning speed is 5°/min. 
Finally, the iNano nanoindenter (KLA-Tencor, The United States) with a load of 50mN was used to 
obtain the mechanical properties of the samples along the upper, middle, and bottom of the section. 
2.3. Calculation 
Molecular dynamics (MD) simulation: An radius 15 A iron sphere with an initial velocity apply an 
initial velocity (Vo) to the LSMP based on the following equations:[17, 18] 
P?(r)=ml,a/2t(a +1) (1) 
Yj=2P(r)t/m (2) 
where J is the nominal laser intensity, t is the pulse duration, m is the mass of iron sphere, P is the 
pressure generated during laser irradiation, and a is the factor to account for the increase of plasma 
thermal energy by the absorbing laser irradiation. A typical velocity of 1 kms | is applied, which yields 
a peak pressure of 0.33Gpa based on the formula (10) in the present case. CoCrFeNi HEAs y and iron 
sphere are single crystals with FCC and Body-centered cubic (BCC) structures, respectively, their 
lattice constants are 3.56 A and 2.863 A. The modified embedded atom method (MEAM) in C style 
was used in the calculation, Mi et al.[19] calculated the potential function suitable for HEAs, which 
was used in this simulation. And the NVT ensemble was applied, warming up to 2000 K with a 0.001 
ps time step. 
Based on the finite element idea, a three-dimensional numerical model of the interaction between 
the laser and the melt pool is established to predict the temperature field during SLM process. In order 
to not affect the calculation quality, the following assumptions are put forward: 1. The plasma and 


metal vapor are ideal gases; 2. The effect of the shielding gas on the welding process is ignored. 3. 


Computational fluids are laminar and Newtonian incompressible fluids. 


Conservation equations: The fluid flow and heat transfer process in the LAM process mainly 


include the conservation equations of mass, momentum, and energy: [20, 21]. 


Conservation of mass: 


O = 
= (p)+V(pV) =0 (3) 
Ot 
Conservation of momentum : 
ð = eoe = = 
a TVN) = MUN) =P= E TAE (4) 
Conservation of energy : 
0 = 
P» (ph) +V(pVh) = V(KVT) (5) 


where p is the fluid density, V is the fluid velocity, t is the time; u is the fluid viscosity, P is the pressure, 
K is the isotropic permeability, g is the gravity; T is the temperature, h is the material enthalpy, and k is 
the thermal conductivity. 


Laser heat source is seted as Gauss volume heat source model:[22] 


(6) 


Quaser = hays Q ex ae. 2 y) | 


mR, H(1-e”) R; log(H / 2) 
where Q is laser power, H is the height of the laser heat source, a,,, is the absorption rate of the 
material to the laser heat source, and Rọ is the effective radius of the laser beam. In theory, more than 
95% of the total energy of the laser beam is concentrated in the area with Roas the radius. 


Driving force: For recoil pressure, a widely accepted model is expressed by:[23] 


P. =0.54P, exp| AH, da (7) 
RTT, 


where Pois the atmospheric pressure, R is the universal gas constant, 4H, is the vaporization enthalpy 


of the liquid metal, T,y is the liquid-gas equilibrium temperature of the material under a certain 


pressure, A Vzy is the difference between the specific volumes of liquid and gas, and Vzy is the specific 


volume of the gas. The surface tension can be expressed as a function of temperature:[20] 


dy 
=y, +— (T-T, 8 
Y trl a) (8) 


Marangoni force formed on the surface of fluid which can be expressed by:[20] 


__ y oT 


G=- (9) 
l 

í OT Ox, 

where T is the liquid surface temperature, Tọ is the reference temperature, and yọ is the surface tension 
at Tọ temperature, Fy-i is the Marangoni force formed by the surface tension gradient in the i direction, 
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y is the surface tension, anı indicates that the direction of the force is opposite to the gradient 
direction of the surface tension. 
Pulse pressure: The pulse pressure is the force formed by accelerating the impact on the surface of 


the molten pool in the vertical direction, and its direction is vertically downward. According to 


theoretical analysis and experimental verification, the pulse pressure formula is as follows:[24] 


a 
P.,, (GPa) = 0.01,| Z JI (10) 
mx (GPa) PEA 


where Pmax is the pulse pressure, a is the efficiency of the interaction, Z is the reduced shock impedance, 


ly is the laser power density. 


Geometric model and boundary conditions: As shown in Fig.2, the computational domain is 30 


mm x 7 mm x 10 mm with the grid resolution of 0.15 mm, 5 mm fluid region represents the HEAs 


powder, and 5 mm void region represents the atmospheric environment and serves to observe the free 


surface. For each boundary, the energy and pressure boundary balance conditions are required. The 


implicit SOR and VOF algorithms are used to obtain and update free surfaces through self-defined 


functions, and boundary conditions are imposed on all meshes until the end of the calculation. 
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Fig.2. The schematic diagram of computational domain. 


3. Results and discussion 


3.1. Surface topography 


To explore the evolution of LSMP, real-time images of the melt pool under the two processes were 


collected as shown in Fig. 3 and Video 1 (Supplementary material). Fig. 3(a) records the real-time 


process of the melt pool in the SLM process. CW laser input fully melts the CoCrFeNi powder and 


gradually converges into a plump melt pool in the front of the solidification zone. In the LSMP process, 


the melt pool exhibits different characteristics. In Fig. 3(b), the plasma clusters induced by the pulsed 


laser can be observed. The most important thing is that the pulsed laser drives and changes convection 


in the melt pool through the shock wave. In the transient process of 3.87ms, the surface area of the melt 


pool expands firstly after each laser shock, 
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Fig.3. Real-time images of melt pool evolution: (a) SLM treatment; (b) LSMP treatment. 


which is beneficial for the escape of heat in the melt pool and increases the cooling rate. Then, the 
surface tension makes the melt pool return to its original state to wait for the next pulsed laser. The 


melt pool cools faster throughout the manufacturing process. 
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3D profile data also confirmed this result, CW laser input fully melts the CoCrFeNi powder and 


results in an arc-shaped melt pool as shown in Fig. 4(b), but this will cause the heat in the melt pool to 


fail escape quickly. Combined with Fig. 4(c), the most striking feature of the LSMP process is that the 


profile of the melt pool has changed, and the size of the melt pool has become wide and flat. The height 


of the melt pool has decreased and the width of the bottom has increased, but this has not damaged the 


molding process. On the contrary, this phenomenon indirectly indicates that the cooling rate of the 


being shaken melt pool is increased, and it has solidified before returning to its original arc-shaped 


state. Owing to the conventional cooling method the heat source exists in the center of the chill layer is 


broken in the LAM process. 


(a) 8. 400mm 


Z 
100 Laser shock L 


v v 


1190 pm 


1323 um 


Height of the melt pool (um) 
Height of the melt pool (um 


L 1 1 L 1 L 


: 3 -200 
0 200 400 600 800 1000 1200 200 400 600 800 1000 1200 1400 
Width of the melt pool (um) Width of the melt pool (um) 


Fig.4. 3D topography and section profile of melt pool under different processes: (b) SLM treatment; (c) LSMP treatment. 


3.2. Evolution of cooling rate 
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Changes in melt pool morphology are attributed to cooling rate. As a result, The CoCrFeNi HEAs 
model was established in MD as shown in Fig. 5(a). Fig. 5(c) shows the liquid melt pool oscillates after 
being stressed, obvious deformation occurs at the top and spreads around. The shock wave propagates 
in three-dimensional directions and the melt pool quickly recovers to its original state. Apparently, 
similar results also appeared. After cooling, the size of the melt pool is different to that of SLM model, 
and becomes wide and flat. The changes in the melt pool profile also agree with the 3D profile data. 
More cooling process under different time steps as shown in Fig. Al (Supplementary material). 
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Fig.5. Melt pool evolution based on Molecular Dynamics: (a) melt pool model; (b) SLM treatment; (c) LSMP treatment. 
The evolution of cooling rate is divided into two parts. From a microscopic perspective, the atoms 


inside the melt pool will accelerate after being forced, which is one of the reasons for the increased 


12 


202203.00067v2 


chinaXiv 


cooling rate. Liu et al.[25] often used the Mean Square Displacement (MSD) to study the motion of 
atoms in a system: 

MSD -25 (t+at)-r (eT (11) 
where ri(t) represents the position of atom i at time t, N is the number of atoms in the system, and 4t 
represents the change in time step. Fig.6 shows The MSD curves of atoms in the melt pool region under 
the two processes. A MSD vector of three quantities involve dx, dy and dz in the LSMP are much larger 
than SLM as shown in Fig.6 (a), atoms inside liquid melt pool oscillates after being stressed and 
spreads around. In addition, there is another feature. The MSD has an obvious jump and a linear 
relationship with time when the pulsed laser acts on the melt pool as shown in Fig.6 (b), indicating the 
movement ability of the atoms is particularly strong at this time, which will accelerate the heat 
dissipation. As the temperature decreases, the atoms in the melt pool area begin to arrange in a regular 
and orderly manner. When the solidification temperature is reached, the MSD curves tend to be 
horizontal, because they have formed a regular arrangement, only a slight relaxation at their 
equilibrium position. More results of MSD under different shock cases as shown in Fig. A2. 


(Supplementary material). 
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Fig.6. The MSD of atoms in the melt pool, (a) the squared dx,dy,dz directions displacements, (b) the total squared displacement. 


From a macro perspective, the melt pool is driven and accelerated to flow, which is another reason 


for the increased cooling rate. The internal melt pool model is built to understand as shown in Fig.7, the 


CW laser energy input forms a beam of force, which moved down to the bottom along the vertical 


centerline of the melt pool, and a symmetrical convection ring appears inside the plump liquid melt 


pool. Then it surges from both sides of the melt pool to the surface, forming the typical Marangoni 


effect. Khairrallah et al.[26] demonstrates the significant effect of the recoil pressure and Marangoni 


convection in laser powder bed fusion (L-PBF). However, the heat source stays in the middle and lower 


region and cannot escape as shown in Fig.7(a), which is the reason why high thermal stress is always 


present. Happily, this problem has been solved due to the assistance of pulsed lasers. It can be observed 


from Fig.7(b) the depression appears in the top region because of the applied force, the original 


symmetrical convective ring is changed, and forming multiple convective loops at the top and bottom 


regions, respectively, which accelerates the flow, and prompting the heat source to move and dissipate 


quickly. 
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Fig.7. Computed temperature and vectors distribution under different processes: (a) SLM treatment; (b) LSMP treatment. 
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The temperature-time relationship of the internal melt pool model is extracted as shown in Fig.8. 
Clearly, different temperature profiles appeared under the two processes. The cooling rate C=4T/S" per 
unit time is calculated, AT is the temperature change, S” is the unit time 1s. After reaching the 
temperature peak, LSMP solidifies rapidly with a cooling rate of -4.84x10°K/s"', which is better than 


SLM. Therefore, LSMP is beneficial to increase the cooling rate. 


1800 


—— SLM ——LSMP 
C=AT/S'=-4.34*107 


C=AT/S'=-4.84x107 


Period (s) 


Fig.8. The computed curve of melt pool between temperature and time. 


3.3. Microstructure analysis 

Based on the layer-by-layer manufacturing principle of LAM, CoCrFeNi powder is remelted and 
solidified due to heating, resulting in material structure changes, and the growth morphology of grains 
is closely related to the cooling rate. Under the SLM process, the cross-section of the melt pool shows a 
non-equilibrium structure composed of columnar dendrites and equiaxed dendrites as shown in Fig.9(a). 
From the local bottom boundary enlarged view of the melt pool, there is a good bond between the HEA 
and the substrate, indicating the formed sample has a good tensile ability. Columnar dendrites were 
observed at the bottom of the melt pool, growing upward perpendicular to the substrate boundaries. In 
the LAM process, the laser irradiation area is the center of the heat source, and the processed area 


belongs to the cooling layer, then a high-temperature gradient is formed between them, which driven 
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the columnar dendrites to grow upward along the temperature gradient from the bottom of the melt 


pool as shown in Fig.9(b). Further, the solid and liquid lines inside the melt pool are promoted to show 


a directional solidification trend from the inside to the outside and the bottom to the top. When 


reaching the middle and upper region of the melt pool, a lower temperature gradient due to the distance 


from the heat source cannot drive the epitaxial growth of columnar dendrites. At the same time, the 


EDS results of the columnar dendrites in Fig.9(b) area show that Cr, Fe, Co, and Ni elements are 


uniformly distributed in the sample as shown in Fig.9(c), and no segregation defects appear. Lin et 


al.[27] also found columnar dendrites are one of the typical features of LAMed HEAs. 
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Fig.9. Samples under SLM treatment: (a) The microstructure of the melt pool section, (b) The local SEM image of columnar 
crystals, (c) The EDS result of in (b); Samples under LSMP treatment: (d) The microstructure of the melt pool section, (e) The 


local SEM image of equiaxed crystal, (f) The EDS result of points A and B in (e). 


In Fig.9(d) and (e), it is found the melt pool is an equilibrium structure composed of equiaxed 


dendrites. Apparently, a Columnar-Equiaxed Transition (CET) occurs. Hunt[28] proposed a model that 


when the volume fraction of equiaxed grains is ®>0.49, complete equiaxed dendrite growth occurs, 


and when the volume fraction of equiaxed grains is ®<0.0066, the structure is considered to be 


completely columnar dendrites. More interestingly, Gaumann et al.,[29] proposed a revised model: 
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Cana a ay 
where G is the temperature gradient, V is the solidification rate, nucleation number NO=2x10'°/m’, 
equiaxed grain volume fraction ®, a and n are constants related to materials, n=3.4 and a 
=1.25x10°(K*4/m s), set ® to a lower critical value ®c=0.0066, Gan et al.[30]also obtained a 
judgment based on the G’/V ratio, that is, when the inside of the melt pool meets the following 


conditions, completely columnar dendrites occur: 


n 
Gar (13) 
V 


where K=2.7x10"(K**/m*4 s), this microstructure relationship is crucial between the limited 
solidification speed and temperature gradient as shown in Fig.10 (a). Obviously, LSMP has changed the 
directional solidification mode under the same G condition. As Fig.10 (b) shows the liquid melt pool 
oscillates from top to bottom after being stressed, and spreads out to form multiple convective loops 
inside melt pool. Besides, the expansion of the surface area leads to an increased cooling rate together 
as shown in Fig.8, which directly accelerated the solidification rate V. As a result, the nucleation rate of 
the grains increases, and the growth rate decreases. Therefore, G'/V<K occurs in the LSMP process, 
the recrystallization mode changes to completely equiaxed dendrite growth, replacing the original 
mode. In Fig. 9(f), the dendrite point A and the interdendritic point B exhibit equal elemental 
compositions of Cr, Fe, Co, and Ni. Abolkassem et al.[31] shows the formation of Cr phase clusters and 
Cr,O; will lead to the reduction of Cr elements during the forming process, while other Cr, Co, and Ni 


elements are likely to form solid solutions with Fe, so the Fe element composition will be high. 
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Fig.10. (a) Microstructure selection map for solidification morphology as a function of G and V,[29] (b) Columnar-Equiaxed 


Transition (CET) occurs in LSMP treatment. 


3.4. Phase composition 


Figure 11 shows the XRD patterns of CrFeCoNi powder and CrFeCoNi HEAs under different 


processing. The three typical peaks of the FCC structure, namely (111), (200), and (220), dominate in 


the range from 20° to 90°, and the results show a stable FCC structure is achieved through all 


processing conditions. Besides, Ni, Cr also forms a BCC solid solution with Fe, Abolkassem et al.[31] 


also obtained the same results. Under the LSMP process, the (111) peak shifts slightly to high angles as 


shown in Fig.11(a), which is caused by the decrease in interplanar spacing. In the orthorhombic crystal 


system, the crystal plane spacing dj, is determined by the crystal plane index (h k 1), and the formula 


for calculating the crystal plane spacing is: 


Ii 


For the cubic crystal system, since a=b=c, the above formula can be simplified to: 


a 


VR? +k? +E 


19 


du= 


(14) 


(15) 
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After peak matching, it is obtained that the lattice constants a of the HEA under the two conditions are 
3.55 A and 2.87 A, respectively, and the interplanar spacing corresponding to the (111) crystal plane is 
calculated from the formula (15) to be 2.05 A, respectively. and 1.66A, as shown in Fig.11(b), at the 
macroscopic stress level in multiple crystal scales, the compressive stress brought by LSMP leads to a 
smaller interplanar spacing, which is shifted to a higher angle in the XRD pattern. Gonzalez et al.[32] 
shows the appearance of residual compressive stress can also effectively inhibit the propagation rate of 


hot cracks, which is beneficial. 
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Fig.11. (a) X-Ray Diffraction patterns and (b) interplanar spacing and lattice constant of samples under different treatments. 


3.5. Nano-indentation analysis 


Changes in microstructure will directly affect the mechanical properties of materials. Under the 


same experimental loading conditions, the nanoindentation results show different positions inside the 


melt pool also exhibit different displacements. In Fig.12 (a) and (b), the displacements of the upper and 


middle region of the melt pool after LSMP treatment are 879 nm and 909 nm, respectively, which are 


better than those of SLM, indicating the deformation resistance of the material is improved. Because it 


is located in the strong influence area of the laser shock wave, the grain structure inside the melt pool 


under LSMP is an equiaxed dendrite, while those of SLM are columnar dendrites. Under the same load, 


the isotropic characteristics of equiaxed dendrites are significantly stronger than that of columnar 
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dendrites in terms of deformation resistance. At the bottom of the melt pool, the opposite trend occurs 
due to the distance from the laser shock wave as shown in Fig.12(c). Since the columnar dendrites are 
tightly bound to the substrate, they will exhibit better combined tensile strength when hard objects are 
pressed in. However, the results of nano hardness show the upper, middle, and bottom of the melt pool 
after LSMP treatment are 3.09GPa, 2.83GPa, and 4.67GPa, respectively, as shown in Fig.12 (d), which 
are all better than those of SLM, indicating the CET under LSMP brings beneficial improvements to 


material properties. 
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Fig.12. Nanoindentation load-displacement in different areas: (a) Top, (b) Middle, and (c) Bottom area of the melt pool; (d) 


Hardness in different areas. 


4. Conclusions 

In the process of LAMed high-entropy alloys, the online regulation of the melt pool is particularly 
important. Anew LSMP method moves the LSP from the solid surface to the liquid melt pool is studied. 
Combined with numerical simulation to understand the action mechanism between laser shock and 
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melt pool, the new process is systematically verified from the three perspectives of forming process, 
microstructures, and mechanical properties. The results show the necessity of introducing laser shock, 
and some important conclusions are summarized as follows: 

1. Revealing the convection mechanism of the melt pool. The plasma induced by the pulsed laser 
acts on the melt pool, which promoted the transformation of melt pool flow mode from Marangoni 
effect to multi-convective loop, and successfully increased cooling rate by accelerating melt pool flow. 

2. The microstructure evolution law of CrFeCoNi HEAs was obtained. After the action of the 
pulsed laser, the solidification rate of the melt pool was increased, the microstructure undergoes CET 
transformation, and macroscopically, the nano hardness of the top, middle and bottom regions is 
improved, and no segregation defects occur. 

3. Realizing non-destructive real-time fabrication of CrFeCoNi HEAs with stable FCC structures, 
and the reduction of (111) interplanar spacing proves that laser shock introduces beneficial compressive 
stress. 
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Appendix A. Supplementary material 

This document provides supplementary information to “Study on laser shocking of melt pool in 
Laser Additive Manufacturing of FeCoCrNi High-Entropy Alloys” 
1. Cooling process 

The Molecular Dynamics (MD) simulation of the Laser Shocking of Melt Pool (LSMP) goes 
through five stages: relax, heating (2000K), constant temperature, shocking and cooling. Before 
entering the cooling stage, the shocking module is removed to observe the liquid surface oscillation 
process. The total step of the cooling is 40ps as shown in Fig.Al, and the SLM process has no shock 
module. Clearly, the liquid melt pool oscillates after being stressed, the shock wave propagates in 
three-dimensional directions as shown in Fig.A1(b). At 4ps, the oscillated molten pool diffused to the 
solid phase and the substrate respectively. Until 20ps, the temperature of the melt pool decreased and 
began to gradually solidify. After cooling, the shape of SLM model has no significant changes. But , 
the shape of LSMP model does not recovers to similar state and becomes flat. The changes in the melt 
pool profile also agree with the 3D profile data. This phenomenon indirectly indicates that the cooling 
rate of the being shaken melt pool is increased, and it has solidified before returning to similar state of 


SLM model . 
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Fig.A1. The cooling process of melt pool based on Molecular Dynamics: (a) SLM treatment; (c) LSMP treatment. 


2. Different shock cases. 


More shock cases including 1.5, 2 and 2.5 km/s were simulated as shown in Fig.A2 (a). The results 


show there is not more benefits when the larger the shock wave acts on the melt pool. Too small will 


not work, too large will destroy the melt pool. In our previous experiments, the laser energy 3J will 


produce beneficial oscillations effect, that is, the case of 1km in MD. At this time, the MSD of the 


atoms in the melt pool is the largest, indicating that the atomic movement ability is particularly strong, 


which will accelerate the escape of the heating source. However, others MSD curve is obviously less 


24 


202203.00067v2 


chinaXiv 


than the case of Ikm. A MSD vector of three quantities involve dx, dy and dz are also displayed 


respectively as shown in Fig.A2. 
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Fig.A2. The MSD of atoms in the melt pool, (a) total squared displacement under different shock conditions; the squared 
dx,dy,dz directions displacements: (b) 1.5 km/s, (c) 2 km/s, and (d) 2.5 km/s. 
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